NASA/TM— 2001-211098 



Friction Factor Characterization for High- 
Porosity Random Fiber Regenerators 

Lanny G. Thieme 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 


August 2001 


Acknowledgments 


The author would like to express his appreciation for the support provided by Andrea Clary and the NASA Glenn 
Flow Calibration Laboratory, operated by Cortez III Service Corporation (presently InDyne, Inc.), which completed 
the flow testing on the regenerator samples. Also, the author would like to thank Maury White and Songgang Qiu 
of Stirling Technology Company (STC) and Steve Welty, formerly of STC, for their efforts in providing the 
regenerator samples, the flow test fixture, and the desired flow ranges. This work was performed in 
part under Interagency Agreement SAA3-198 with the Department of Energy, between 

November 1997 to March 1999. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / / gltrs.grc.nasa.gov/GLTRS 



FRICTION FACTOR CHARACTERIZATION FOR HIGH-POROSITY 
RANDOM FIBER REGENER ATORS 


Lanny G. Thieme 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland. Ohio 44135 


Summary 

NASA Glenn Research Center, the Department of Energy (DOE), and Stirling Technology 
Company (STC) of Kennewick, WA are developing a Stirling convertor for a high-efficiency 
Stirling Radioisotope Power System to provide spacecraft on-board electric power tor NASA 
deep space missions. Stirling is also now being considered for unmanned Mars rovers, especially 
for missions of long duration. STC is developing the 55-We Technology Demonstration 
Convertor (TDC) under contract to DOE. NASA Glenn is conducting an in-house technology 
project to assist in developing the convertor for readiness for space qualification and mission 
implementation. NASA Glenn is also evaluating key technology issues through the use of NASA 
Phase II Small Business Innovation Research (SBTR) contracts with STC and providing technical 
consulting for the DOE effort under an Interagency Agreement. 

As part of the Interagency Agreement, NASA Glenn performed a review of the TDC design. 

The HFAST computer code was used to compare projected convertor performance to that 
predicted by the GL1MPS code used by STC to design the convertor. These comparisons 
showed a large difference in the regenerator pressure drop losses due to significantly different 
regenerator friction factors used by the two codes. Neither code s friction factor correlations 
appeared to be based on regenerator test data at the regenerator porosities that were being 
considered for the design (90-96%). Therefore, NASA Glenn recommended that steady-flow 
pressure drop tests be conducted to better determine the friction factor for these high-porosity 
regenerators. 

STC fabricated a flow test fixture and three random fiber regenerator test samples, one each at 
approximately 80%, 88%, and 96% porosities. The flow tests were then completed by the 
NASA Glenn Flow Calibration Laboratory, and the data reduced to Reynolds number and 
friction factor. The results showed that the 80% and 88% porosity samples had similar 
characteristics while the 96% porosity sample had significantly higher friction factors for given 
Reynolds numbers compared to the samples with lower porosities. Comparisons were also made 
between the test data and existing correlations. STC used this data to derive a modified 
regenerator friction factor correlation for use in GLIMPS for porosities greater than 88 /o. Using 
this new correlation, the final optimized regenerator design porosity was reduced from 96% to 
90%. 
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It is important to note that no adjustments were made to the heat transfer correlation used in 
GUMPS for the higher regenerator porosities. Heat transfer test data would be very worthwhile 
to obtain but are difficult to measure. It was felt that the straightforward steady-flow pressure 
drop tests would give at least some idea of the similarity of high-porosity regenerators to those 
for which both friction factor and heat transfer data exist. In the future, an existing NASA Glenn 
regenerator test rig may be used to perform heat transfer testing on high-porosity random fiber 
regenerator samples. 

Introduction 


NASA Glenn Research Center, the Department of Energy (DOE), and Stirling Technology 
Company (STC) of Kennewick, WA are developing a Stirling convertor for an advanced 
radioisotope power system to provide spacecraft on-board electric power for NASA deep space 
missions. Stirling is being evaluated as an alternative to replace RTG’s (Radioisotope 
Thermoelectric Generators) with a high-efficiency power source and has been identified for 
potential use on the Europa Orbiter and Solar Probe missions. Stirling is also now being 
considered for unmanned Mars rovers, especially for missions of long duration. The efficiency of 
the Stirling system, in excess of 20%, will reduce the required amount of isotope by a factor of at 
least 3 compared to RTG’s. This significantly reduces radioisotope cost, radiological inventory, 
and system cost and provides efficient use of scarce radioisotope resources. 

STC has designed, fabricated, and completed first testing of the 55-We Technology 
Demonstration Convertor (TDC) under contract to DOE [1,2]. TDC’s are now being tested 
singly and in dynamically-balanced opposed pairs by STC. The design power of 55 We and an 
efficiency of about 25% have been demonstrated for TDC operating conditions of 650°C hot-end 
temperature and 120°C cold-end temperature. The TDC's have been mapped over a range of 
temperatures and strokes. NASA Glenn is providing technical consulting for this effort under an 
Interagency Agreement with DOE. 

NASA Glenn is conducting an in-house technology project to assist in developing the convertor 
for readiness for space qualification and mission implementation [3]. As part of this effort, the 
TDC has been characterized for electromagnetic interference/electromagnetic compatibility 
(EMI/EMC) and also successfully passed launch environment random vibration testing at 
qualification levels while operating at full stroke and full power. An independent performance 
verification and mapping of multiple TDC’s is now being done. Heater head life assessment and 
NdFeB magnet aging characterization tasks are underway. Substitute organic materials for the 
linear alternator for use in a high radiation environment have been identified and have been 
incorporated by STC in TDC’s recently built for GRC. Electromagnetic and thermal FEA for the 
linear alternator are also being conducted. 

As part of the overall radioisotope Stirling development, NASA Glenn is evaluating key 
technology issues through the use of two NASA Phase II SBIR contracts with STC [4], Under 
the first SBIR, STC demonstrated a synchronous connection of two thermodynamically 
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independent Stirling convertors. Synchronization was achieved over a wide operating range, and 
a 40 to 50 fold reduction in vibrations compared to an unbalanced convertor was shown. System 
operation was also demonstrated with the synchronized convertors feeding a battery charger 
load. This connection method is now' being used to connect the TDC’s. The second SBIR 
contract is for the development of an Adaptive Vibration Reduction System (AV RS) that will 
further reduce vibration levels under normal operating conditions and will add the ability to adapt 
to any changing convertor conditions over the course of a mission. 


As part of the Interagency Agreement. NASA Glenn performed a review of the TDC design. As 
part of this review, the HFAST computer code was used by NASA Glenn to analyze the TDC 
performance. HFAST was developed by Mechanical Technology, Inc. (MTI), in part, during the 
NASA effort to develop Stirling for space power for the SP-100 program. NASA Glenn has 
unrestricted rights to use and distribute HFAST. Comparisons were made to the STC GLIMPS 
computer code predictions for the TDC. GLIMPS was the primary code used by STC for design 
purposes. NASA Glenn also owns a license to use the GLIMPS code and has completed 
previous comparisons between HFAST and GLIMPS. NASA Glenn s past experience in using 
the two codes has not revealed any clear superiority of one code over the other. Although there 
are differences between them, both have been found to give reasonable overall results. 

Overall comparisons of the GLIMPS and HFAST projections for the TDC showed that the two 
codes predicted similar performance with an adjustment for friction factor in the regenerator to 
make the two codes use similar regenerator friction factors. Without this adjustment, HI- AST 
predicted much lower performance than did GLIMPS due to higher pressure drop losses through 
the regenerator. 

The friction factor correlations used in each code were evaluated to understand these differences 
in pressure drop losses. This then led to a recommendation to complete steady-flow pressure 
drop tests to better determine the friction factor for the high-porosity regenerator that is used in 
the TDC design. NASA Glenn performed these tests in the NASA Glenn Flow Calibration 
Laboratory. This report describes the evaluation of the GLIMPS and HFAST friction factor 
correlations, the steady-flow testing and data reduction procedures, and the test results and 
comparisons to existing correlations. 

Comparison of the HFAST and GLIMPS Friction Facto r Correlations 

STC was considering porosities of about 90-96% for the TDC random fiber regenerator design. 
Both GLIMPS and HFAST include a porosity dependence in their friction factor correlations. 
However, as shown by figures 1 and 2, HFAST has a much higher sensitivity to porosity than 
does GLIMPS, especially for screen regenerators. The agreement between the two codes is 
better for 94% random fiber regenerators as shown in figure 3, but the codes still vary by about 
40%. The friction factor correlations for both codes agree well at 78% porosity for both types of 
regenerators. 
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Figure 1 - GLIMPS and HFAST Friction Factor (f) vs. Reynolds Number (Re) for 
94% Porosity Screen Regenerators. 
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Figure 2 - GLIMPS and HFAST Friction Factor (f) vs. Reynolds Number (Re) for 
78% Porosity Screen Regenerators. 


Gedeon, in his final report on regenerator pressure drop and heat transfer testing in oscillating 
flow done at Ohio University [5], concluded that he could not determine any porosity 
dependence for the range of screens and felts (random fiber) that were tested in that project. 
These tests covered porosities of 62-78% for screens and 69-84% for felts. The HFAST manual 
states that its screen regenerator friction factor is based on Kays and London data and appears to 
be the same as used in an earlier version of GLIMPS. Kays and London show that their test 
data were taken for a range of porosities of 60-83% for screens [6]. GLIMPS (version 4.0) uses 
relationships for screens and random fibers that are based, in part, on earlier testing at Ohio 
University [7], Pressure drop tests for steady and oscillating flow were also completed by 
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Figure 3 - GLIMPS and HFAST Friction Factor (f) vs. Reynolds Number (Re) for 
94% Porosity Random Fiber Regenerators. 


Sunpower, Inc. with the same rigs that were later moved to Ohio University [8]. These tests 
included regenerator porosities of 61-68% for screens and 80-84% for random fibers. 

So it appears that neither code's friction factor correlations are based on regenerator test data at 
the regenerator porosities that were being considered for the design (90-96%). Due to this 
concern, NASA Glenn recommended that STC consider running a steady-flow pressure drop test 
on a high-porosity regenerator sample to get good empirical data on the friction factor for these 
regenerators. Gedeon indicated that steady-flow testing should be adequate for determining 
friction factors for typical regenerators [7]. Heat transfer test data would also be very important 
but are much more difficult to obtain. It was felt that the straightforward pressure drop tests 
would give at least some idea of the similarity of high-porosity regenerators to those for which 
both friction factor and heat transfer data exist. 


Steadv-Flow Pressure Drop Tests 

STC agreed with the recommendation to conduct steady-flow pressure drop tests to determine 
the friction factor for a high-porosity regenerator sample. It was determined that the Flow 
Calibration Laboratory at NASA Glenn could do these tests in a quick, economical, and accurate 
manner. 

STC selected and fabricated three random fiber regenerator test samples, one each at 
approximately 80%, 88%, and 96% porosities. The 80% porosity sample was chosen to give a 
comparison to the correlations used in GLIMPS and HFAST at a porosity level where empirical 
test data exist. The wire diameter used in all samples was 0.0009 in (22 microns). Each test 
sample consisted of two disks that were stacked in the sample holder. STC provided the length, 
diameter, and porosity of each disk. The values of diameter and porosity for the two disks in 
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each sample were averaged and the lengths of the two disks were added to arrive at the following 
sample characteristics: 

Sample Ave. diameter, cm (in) Ave. porosity Length, cm (in) 

80%, disks 1 ,2 2.22 (0.875) 0.797 2.57(1.011) 

88%, disks 3,4 2.22 (0.873) 0.878 2.53 (0.995) 

96%, disks 5.6 2.21(0.872) 0.957 2.42 (0.952) 

STC also designed a test fixture, with inputs from NASA Glenn, and then fabricated this fixture 
to hold the samples for testing. The test fixture and the three test samples are shown in figure 4. 
The pressure taps can be seen on each half of the test fixture. A coarse wire screen was used on 
each halt of the test fixture to help straighten the flow on each end of the regenerator sample. 
Photomicrographs of sections of each of the three test samples are shown in figure 5. 

Consideration was given for flow testing with helium, nitrogen, or air. The TDC uses helium at 
about 2.41 MPa (350 psia) absolute mean pressure; the maximum flow rate in the regenerator 
predicted by GL1MPS was about 3.6 g/sec (0.008 lbm/sec). Based on discussions with the 
NASA Glenn Flow Calibration Laboratory and STC, it was decided that testing with air at about 
0.69 MPa (100 psia) absolute inlet pressure provided the best combination of meeting the desired 
flow rates and producing the data in a timely, economical manner. STC provided their desired 
flow rates for testing. These were chosen to both cover the range of Reynolds numbers as 
appropriate for the TDC and to include higher flow rates that may be applicable to future 
convertor designs. 



Figure 4 - Regenerator test samples and flow test fixture. 
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Figure 5 - Random Fiber Regenerator Samples at 50x Magnification. 

The flow tests were then completed by the NASA Glenn Flow Calibration Laboratory. The test 
data are included in Appendix A. A test was also run with no sample to determine the effect of 
the coarse screens that were used as flow straighteners. There was no measurable pressure drop 
for this test until the flow rate reached 5.4 g/sec (0.012 lbm/sec). At the highest flow rate, 

1 1 .8 g/sec (0.026 lbm/sec), the pressure drop was only 0. 1 7 kPa (0.025 psid). It was felt that 
this pressure drop was negligible, and it was ignored in the data reduction. 

Figure 6 shows the measured steady-flow pressure drop as a function of mass flow rate for each 
of the three samples. Testing was done with air at 0.69 MPa (100 psia) absolute inlet pressure 
for all samples. Flow conditions were chosen to include the range of Reynolds numbers that 
were expected to occur in the actual convertor. Tests were also run in each of the two different 
flow directions for the 96% porosity sample. These results are shown in figure 7 and can be seen 
to be essentially identical. Flow directions 5-6 and 6-5 refer to the regenerator disk numbers 5 
and 6, which were the two 96% porosity disks that were used in this sample. 
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Mass Flow, g/sec 

Figure 6 - Pressure Drop vs. Mass Flow for Random Fiber Regenerator Samples. 



Mass Flow, g/sec 

Figure 7 - Pressure Drop vs. Mass Flow for 96% Porosity Random Fiber Regenerator 
Sample for Different Inlet Pressures and Different Flow Directions. 


A test was also run at 0.38 MPa (55 psia) absolute inlet pressure for the 96% porosity sample to 
check the effect of AP/P, and these pressure drop results are compared to those for the 0.69 MPa 
(100 psia) absolute inlet pressure in figure 7. An expansion factor correction based on the AP/P 
[9], calculated using inlet pressure, was made on one of these sets of data to correct it back to the 
data at the other pressure level. The thought was to then use this type of correction on any data 
points that had high ratios of AP/P. However, it was found that similar results were obtained by 
using the average of the inlet and outlet pressures when calculating the density of the fluid so this 
method was used for all data and assumed to be satisfactory for any data points that had a high 
ratio of AP/P. In general, it was preferred to maintain the AP/P below 0. 1 . Appendix A shows 
that most of the test data had a AP/P of less than 0.06. 
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The measured flow data were provided to STC. STC and NASA Glenn first independently 
reduced this data to friction factor versus Reynolds number and then resolved any significant 
differences. The Reynolds numbers and friction factors as calculated by NASA Glenn are 
included in Appendix A. These values were obtained from the test data by using the following 
equations: 

Re = mdot * D h /(p * A samp | C * p) 


D h = p * d w / ( 1 -p) 


where: 


Re - Reynolds number 

mdot - mass flow rate, lbm/sec 

D h - regenerator sample hydraulic diameter, ft 

p - fluid absolute viscosity, lbm/ft-sec 

Asampic - total cross-sectional area of the regenerator sample, ft : 

p - regenerator sample porosity 

d w - regenerator sample wire diameter, ft 

f = AP * D h * (Ajampic * p) 2 * P * 2g c * 144 / (L * mdot 2 ) 

where: 

f - friction factor 
AP - pressure drop, psi 

Dh - regenerator sample hydraulic diameter, ft (as defined in previous 
equation) 

Asampic - total cross-sectional area of the regenerator sample, ft 2 
p - regenerator sample porosity 
p - fluid density, lbm/ft 3 

gc - gravitational conversion factor = 32.2 lbm-ft/(lbf-sec“) 

L - regenerator sample length, ft 
mdot - mass flow rate, lbm/sec 

Friction factor versus Reynolds number (based on NASA Glenn’s data reduction) is shown for 
each of the samples in figure 8. The curves for the 80% and 88% porosity samples are very 
similar while the friction factor is significantly higher for the 96% porosity sample. Figure 9 
shows the friction factor versus Reynolds number results for the two different pressure levels 
run for the 96% porosity sample. They are identical for most of the Reynolds number range 
with significant variations at the lower Reynolds numbers. These variations appear to be caused 
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‘88% Porosity 
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Re 

Figure 8 - Friction Factor (f) vs. Reynolds Number (Re) for Random Fiber 
Regenerator Samples. 



Figure 9 - Friction Factor (f) vs. Reynolds Number (Re) for 96% 
Regenerator Sample for Different Inlet Pressures. 


-A- 0.38 MPa 
“•“0.69 MPa 


Porosity Random Fiber 


by errors in the measurement of pressure drop as the values of the pressure drop are very low for 
these data points. 


Estimated errors for all data are shown in figure 10. Errors were calculated based on the following 

[ 10 ]: 

W R = [((3R/5X|)W|) 2 + ((dR/dx 2 )W 2 ) 2 + ... + ((dR/dx^Wn) 2 ] 1 ' 
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♦“80% Porosity 
■A— 88% Porosity 

96% Porosity - 0.69 MPa 
96% Porosity - 0.38 MPa 


Figure 10 - Percent Error in the Friction Factor (f) vs. Reynolds Number (Re) for Random Fiber 

Regenerator Flow Tests (the error in Reynolds number is 1.0% of reading for all points). 


where: 


w R - uncertainty of the calculation, units of the calculation 

dR/3x,, 3R/dx 2 , ..., 3R/3x n - partial derivatives of the calculation with respect to each 
independent variable in the calculation 

w i , w 2 , ...» w n — uncertainties of the independent variables, units of the independent 
variables 

The measurement error for the mass flow rate was 1% of reading for all data points. The error in 
Reynolds number was taken as 1% of reading based on just this error in the mass flow rate. 
Uncertainties in the pressure drop, mass flow rate, inlet pressure, and temperature were included 
in the calculation of uncertainty for the friction factor. Figure 10 shows the error in the friction 
factor as a function of Reynolds number for each of the three samples; curves are shown for each 
pressure level for the 96% porosity data. Most of the errors in friction factor are less than 10%. 
Significantly higher errors occur for the low Reynolds number points for the 96% porosity 
curves. These larger errors appear to explain the difference in friction factor for the two different 
pressure levels at these low Reynolds numbers for the 96% porosity sample. 

The friction factor curves for the three samples are compared with three correlations in figures 
1 1-13. The three correlations shown are those for GLIMPS, HFAST, and the latest correlation 
based on test data taken at Ohio University [5]. All three correlations agree well with the test 
data for the 80% porosity sample (figure 1 1 ). Figure 1 2 shows that the GLIMPS correlation and 
the correlation based on test data taken at Ohio University agree well with the test data for the 
88% porosity sample while the HFAST correlation yields a higher friction factor. Finally, figure 
13 shows the same comparison for the 96% porosity sample. Now, HFAST compares the best 
with the test data (data shown is for the 0.69 MPa (100 psia) absolute inlet pressure only) at the 
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Figure 1 1 - Friction Factor (f) vs. Reynolds Number (Re) for 80% Porosity Random Fiber 
Regenerator - Comparison of Test Data with Correlations. 



Re 

Figure 12 - Friction Factor (f) vs. Reynolds Number (Re) for 88% Porosity Random Fiber 
Regenerator - Comparison of Test Data with Correlations. 


low-er Reynolds numbers, 0-300, which was the expected range for the TDC design with a 96% 
porosity regenerator. 


Following analyses of these results, STC decided to use the GLIMPS correlation for all 
regenerator porosities up to 88%. They then derived a curve based on the GLIMPS friction 
factor for 88% porosity and the 96% porosity test data (average of the curves for 0.38 MPa 
(55 psia) and 0.69 MPa (100 psia) absolute inlet pressures) that extrapolates between these for 
porosities from 88-96%. Figure 14 shows the GLIMPS friction factor curve for 88% porosity 
and the two test data curves for the 96% porosity sample over the Reynolds number range of 
0-350. 
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Figure 13 - Friction Factor (f) vs. Reynolds Number (Re) for 96% Porosity Random Fiber 
Regenerator - Comparison of Test Data with Correlations. 



Re 

Figure 14 - Input Curves for STC Correlation for Friction Factor for Random Fiber 
Regenerator Porosities Between 88% and 96%. 


Based on these overall results and recognizing that heat transfer data does not exist for the higher 
porosities, STC decided to limit the design regenerator porosity to a value no greater than 90%. 
STC GLIMPS results with the new friction factor input for 88-96% porosities showed that the 
final design optimized around 90% porosity (with no constraint on the porosity), and this is the 
regenerator porosity that was then used for the TDC. It should be kept in mind that no 
adjustments were made to the heat transfer correlation used in GLIMPS for the higher porosities. 
In the future, an existing NASA Glenn regenerator test rig may be used to perform heat transfer 
testing on high-porosity random fiber regenerator samples. 
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Conclusions 


Prior to the testing discussed in this report, the random fiber regenerator for the DOE/STC TDC 
design was optimizing at porosities of about 96%. Comparisons of the GLIMPS and HFAST 
computer codes showed large differences in the regenerator friction factor correlations for these 
porosities although the correlations compared well for porosities around 80%. It was felt that 
the large differences for the higher porosities were due to the correlations being extrapolated too 
far from existing regenerator test data on which the correlations were based. The highest- 
porosity regenerator flow’ test data from which the GLIMPS and HFAST correlations were 
derived appear to be for porosities of around 84%. Thus, it was recommended to perform flow 
tests on high-porosity random fiber regenerators to better determine the friction factor. 

Steady-flow tests were run for 80%, 88%, and 96% porosity random fiber regenerator samples, 
and the data were reduced to Reynolds number and friction factor. The results showed that the 
80% and 88% porosity samples had similar characteristics while the 96% porosity sample had 
significantly higher friction factors for given Reynolds numbers compared to the lower porosities 
STC used this data to derive a modified regenerator friction factor correlation for use in GLIMPS 
for porosities greater than 88%. Using this new correlation, the final regenerator design 
optimized at a porosity of 90%. 

It is important to note that no adjustments were made to the heat transfer correlation used in 
GLIMPS for the higher regenerator porosities. Heat transfer test data would be very w orthwhile 
to obtain but are difficult to measure. It w'as felt that the straightforward pressure drop tests 
would give at least some idea of the similarity of high-porosity regenerators to those for which 
both friction factor and heat transfer data exist. In the future, an existing NASA Glenn 
regenerator test rig may be used to perform heat transfer testing on high-porosity random fiber 
regenerator samples. 
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Sample: 80 % Porosity - Disks 1,2 

Flow Direction: 1-2 Date: 3 / 3/98 
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Sample: 96% Porosity - Disks 5,6 

Flow Direction: 5-6 Date: 3/3/98 
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100.02 689.6 0.0226 10.3 1.11 7.7 0.011 76.0 24.4 

100.06 689.8 0.0268 12.1 1.47 10.1 0.015 76.0 24.4 

104.52 720.6 0.0164 7.44 0.63 4.3 0.006 76.1 24.5 

101.47 699.6 0.00848 3.85 0.23 1.6 0.002 76.1 24.5 

100.07 690.0 0.00256 1.16 0.05 0.3 0.0004 76.2 24.6 
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